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Fortran  digital-computer  routines  to  calculate  absorption  loss  for  radiowave  propagation  in  the  tropo¬ 
sphere  have  been  developed,  and  machine  plots  of  radar  attenuation  as  a  function  of  range  have  been  made 
for  a  standard  atmosphere,  over  the  frequency  range  100  MHz  to  100  GHz.  The  tropospheric  noise  temper¬ 
ature  has  also  been  calculated  and  plotted.  This  work  is  an  updating  of  earlier  calculations  that  were  less 
accurate,  especially  in  the  frequency  range  above  10  GHz.  Separate  curves  Tor  the  oxygen  and  water-vapor 
components  of  the  absorption  loss  are  provided  in  addition  to  the  curves  for  total  loss;  this  allows  the  loss 
for  values  of  water-vapor  density  other  than  the  standard  7.5  g/m3  to  be  found,  by  applying  a  simple 
multiplicative  factor.  The  computer  routines,  listed  at  the  end  of  the  report,  have  also  been  used  as  part 
of  a  computer  program  to  calculate  the  maximum  range  of  a  radar  system. 
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ABSTRACT 


♦ 


Fortran  digital-computer  routines  to  calculate  absorption  loss  for  radio¬ 
wave  propagation  in  the  troposphere  have  been  developed,  and  machine 
plots  of  radar  attenuation  as'a  function  of  range  have  been  made  for  a 
standard  atmosphere,  over  th?  frequency  range  100  MHz  to  100  GHz.  The 
tropospheric  noise  temperature  has  also  been  calculated  and  plotted.  This 
work  is  an  updating  of  earlier  calculations  that  were  less  accurate,  especially 
in  the  frequency  range  above  10  GHz.  Separate  curves  for  the  oxygen 
and  water-vapor  components  of  the  absorption  loss  are  provided  in  addition 
to  the  curves  for  total  loss;  this  allows  the  loss  for  values  of  water-vapor 
density  other  than  the  standard  7.5  g/m3  to  be  found,  by  applying  a  simple 
muli  iplicative  factor.  The  computer  routines,  listed  at  the  end  of  the 
report,  have  also  been  used  as  part  of  a  computer  program  to  calculate  the 
maximum  range  of  a  radar  system. 
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RADAR/RADIO  TROPOSPHERIC  ABSORPTION  AND 
NOISE  TEMPERATURE 


INTRODUCTION 

The  objective  of  the  work  reported  here  was  to  develop  a  computer  program  for 
calculating  the  absorption  of  radio  waves  by  the  oxygen  and  water  vapor  of  the  earth’s 
lower  atmosphere  —  the  nonionized  region  known  as  the  troposphere.  Resulting  calcula¬ 
tions  have  been  used  to  produce  machine-plotted  curves  of  attenuation  as  a  function  of 
frequency  and  other  parameters.  The  requirement  that  instigated  this  work  was  the  need 
to  know  the  magnitude  of  the  absorption  loss  in  calculating  the  maximum  range  of  a 
radar.  The  computer  subroutines  to  be  described  have  in  fact  been  incorporated  into  a 
larger  program  that  computes  radar  maximum  range  (1).  However,  the  results  are  equally 
applicable  to  radio  systems  of  any  type  that  use  the  troposhere  as  a  propagation  medium. 

The  frequency  range  considered  is  100  MHz  to  100  GHz.  Below  100  MHz  (and 
even  somewhat  above  this  frequency)  tropospheric  absorption  is  negligible.  Above  100  GHz 
strong  water-vapor  resonance  lines  result  in  great  absorption:  this  region  is  not  of  much 
current  interest  for  radar. 

The  work  described  updates  earlier  work  begun  in  1958,  described  in  various  NRL 
reports  and  other  writings  (2-8).  The  original  work  (2,4,6)  resulted  in  curves  plotted 
manually  from  manual  calculations  for  frequencies  from  100  MHz  to  10  GHz.  Later 
(7,8),  the  NRL  NAREC  digital  computer  was  used  to  recalculate  the  absorption  and  noise 
temperature  more  exactly  and  to  extend  the  frequency  range  to  100  GHz.  However  in 
that  work  (as  in  the  earlier  work)  the  “centroid  approximation’’  of  Van  Vleck  (9)  was 
used  to  calculate  the  oxygen  absorption;  hence  the  results  were  not  accurate  i:i  and  near 
the  60-GHz  oxygen  resonances,  that  is,  in  the  region  from  about  50  to  70  GHz.  Also, 
the  so-called  residual  absorption  below  100  GHz  due  to  strong  water-vapor  resonances  at 
higher  frequencies  was  approximated  by  a  formula  suggested  by  Van  V'mok  (10),  but  the 
constant  factor  used  has  subsequently  been  shown  to  give  absorption  values  that  are  too 
low. 


In  the  current  work,  presented  in  tins  report,  the  oxygen  resonances  are  explicitly 
and  individually  included  in  the  calculation,  so  that  the  results  are  valid  in  the  region 
from  50  to  70  GHz.  The  equations  us°d  are  those  of  Meek'  and  Lilley  (11)  as  slightly 
modified  by  Rehcr,  Mitchell,  and  Carte-  (12).  Improved  numerical  integration  formulas 
are  used  for  the  cumulative-absorption  and  noise-tempeature  calculations.  An  improved 
expression  is  used  for  the  residual  effect  of  the  millimeter-wavelength  water -vapor 
resonances,  with  a  constant  based  on  experimental  results  near  100  GHz  reported  by 
Straiton  and  Tolbert  (13).  Also  a  more  sophisticated  equation  for  the  absorption  of  the 
22.235-GHz  water-vapor  line,  due  to  Liebe  (14),  is  used.  Finally  an  improved  model  of 
the  atmospheric  water-vapor  content  is  used.  The  numerical  results  of  these  improved 
calculations  differ  significantly  from  the  previous  results  only  above  about  10  GHz.  Below 
that  frequency  the  differences  are  small. 
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Many  authors  have  made  absorption  calculations,  and  it  is  therefore  pertinent  to  state 
why  it  is  necessary  or  desirable  to  publish  these  of  the  present  report.  Among  others 
the  reports  and  papers  of  Bean  and  Abott  (15),  Hogg  (16),  Schulkin  (17),  Le  Fande  (IK), 
and  Ulaby  and  SUaiton  (19)  may  be  cited,  in  addition  to  those  previously  mentioned. 

This  list  is  not  exhaustive,  in  particular,  it  does  not  include  all  of  the  work  published  in 
other  countries,  some  of  which  is  referenced  by  the  authors  Htcd.  The  published  results 
by  these  authors  have  been  either  limited  to  calculation  of  the  absorption  coefficient 
(absorption  per  unit  distance)  as  a  function  of  frequency  rather  than  calculation  of  the 
cumulated  absorption  along  a  rr.y  path  or  limited  to  calculation  of  the  total  absorption 
for  paths  traversing  the  entire  troposphere.  Some  results  are  limited  to  specific  frequency 
regions  or  absorption  lines.  Moreover  the  curves  presented  by  most  other  authors  cannot 
be  read  with  sufficient  accuracy  to  be  useful  for  practical  calculations  of  radio  or  radar 
system  performance. 

For  radar  detection-range  calculations  the  absorption  between  two  points  within  the 
troposphere  is  usually  needed,  rather  than  total  absorption  through  the  entire  troposphere. 
This  report,  and  its  predecessors,  present  curves  of  the  former  type  (as  well  as  the  latter), 
with  coordinate  grids  and  scales  that  permit  reading  values  to  at  least  two-significant- 
figure  accuracy.  The  curves  of  the  present  report  have  been  plotted  using  the  NRL  Gerber 
plotter  (Gerber  Model  875  automatic  drafting  machine)  and  special  computer  plotting 
subroutines.  The  coordinate  grids  as  well  as  the  curves  are  plotted  in  a  single  operation, 
to  avoid  registration  errors.  The  plotting  accuracy  is  of  the  order  of  0.001  inch  on  the 
original  plots,  which  have  been  photographically  reduced  to  about  1/3  to  1/5  as  large  for 
reproduction  in  this  report. 

To  the  author's  knowledge  results  published  by  others  based  on  computer  calculation 
have  not  been  accompanied  by  listings  of  the  computer  programs  used.  The  programs 
previously  written  for  the  NAREC  computer,  on  which  the  author's  previous  reports 
were  based,  were  written  (by  Maurice  Brinkman,  of  the  NRL  Research  Computation 
Center)  in  the  NELIAC  language,  which  was  never  widely  used  and  is  no  longer  much 
used  at  NRL.  The  newer  routines,  described  in  this  report,  have  been  written  by  the 
author  in  Fortran  for  the  NRL  CDC-3800  computer.  Complete  listings  of  these  routines 
will  be  given  in  the  last  section  of  this  report.  Although  versions  of  the  Fortran  language 
used  with  other  machines  may  difler  slightly,  tne  differences  are  usually  minor,  so  that 
the  routines  can  be  readily  adapted  to  any  computer  that  has  a  Fortran  compiler. 

Another  feature  of  the  curves  presented  in  this  report  is  the  separation  of  the  oxygen 
and  water-vapor  absorption  losses,  in  addition  to  curves  showing  the  combined  absorption 
loss.  The  combined  results  arc  for  a  “standard”  atmospheric  water-vapor  density  of  7.5 
grams  per  cubic  meter.  The  oxygen  component  of  absorption  relatively  stable  and 
insensitive  to  day-to-day  changes  in  the  atmosphere.  For  various  atmospheric  conditions 
however  the  water -vapor  conter.'  may  vary  between  roughly  2  and  20  grams  per  cubic 
meter.  The  water-vapor  absorption  is  (except  tor  a  small  an  !  gcnrrr.il/  i.i>iK  ificant  non¬ 
linearity)  proportional  to  the  water-vapor  density.  Consequently  separate  oxygen  and 
water-vapor  curves  allow  the  water-vapor  component  to  be  adjusted  for  any  water-vapor 
content  other  than  the  standard  7.5  grams  per  cubic  meter  by  simply  multiplying  the 
7.5-gram  absorption  value  by  p/7.5,  where  p  is  the  actual  gram-per-cubic-meter  water-vapor 
content.  This  corrected  water-vapor  absorption  can  then  be  directly  added  to  the  oxygen 
absorption  to  obtain  the  total  absorption. 
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The  curves  presented  for  cumulative  attenuation  along  a  ray  path  are  for  the  rud&r 
case  of  two-way  traversal  of  the  path.  Consequently  the  decibel  attenuations  are  exactly 
twice  as  great  as  they  would  be  for  the  radio-communication  case  of  one-way  propagation. 
The  curves  can  be  used  for  one-way  systems  by  simply  taking  half  of  the  decibel  attenuation 
values  shown  for  the  radar  case. 


ABSORPTION  COEFFICIENT  FOR  OXYGEN 

As  mentioned  in  the  Introduction,  the  absorption  coefficient  for  oxygen  has  been 
calculated  using  equations  given  by  Meeks  and  Lilley  (11),  as  slightly  modified  by  Reber 
et  al.  (12).  The  theory  was  originally  duo  to  Van  Vleck  (9);  later  workers  have  refined 
some  of  the  details.  In  particular  Meeks  and  Lilley  have  formulated  a  refinement  of  the 
dependence  of  the  line-breadth  constant  (which  will  here  be  denoted  Af)  on  altitude 
(pressure  and  temperature),  and  Reber  et  al.  have  made  a  slight  further  refinement  of 
this  factor. 

The  Meeks  and  Lilley  equations  are  a  summation  of  the  absorption  over  many  oxygen 
resonance  lines  in  the  vicinity  of  60  GHz.  These  resonances  occur  for  odd-integral  values 
of  the  rotational  quantum  number  N,  that  is,  for  N  =  1,  3,  5,  ...  .  It  has  been  determined 
(11)  that  the  contribution  of  terms  above  N  =  45  is  negligible;  hence  45  is  the  largest 
value  included  in  the  summation.  For  each  odd  value  of  N  there  are  two  resonance 
frequencies,  denoted  fN+  and  fN_.  These  resonances  are  listed  in  Table  1,  as  given  by 
Meeks  and  Lilley.  The  absorption  at  an  arbitrary  frequency  f  is  the  sum  of  contributions 
from  each  of  these  collision-broadened  resonance  lines,  plus  a  nonrosonani  component. 

The  broadened  resonance  lines  are  assumed  to  have  shapes  defined  by  the  Vnn-Vleek- 
VVeisskopf  formula: 


Fn±  (f)  = 


-  f)2  +  ( Af)2  (fNi  +  f)2  +  (Af)2 


The  nonresonant  contribution  is  of  the  form 


=  Af 
f2  +  (Af)2 


(2) 


The  terms  summed  over  odd  values  of  N'  are  of  the  form 


AN  _  ^N*  +  FN  •  ^N-  +  *0  0  *N  > 


where 


2  N(2N  +  3) 
^N+  N  +  1 


(4) 


2  _  (N  +  1)  (2N  —  1 ) 

M 


(5) 
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Table  1 

Oxygen  Resonance  Frequencies 


N 

fN+  (GHz) 

In—  (GHz) 

1 

56.2648 

118.7505 

3 

58.4466 

62  4863 

5 

59.5910 

60.3061 

7 

60.4348 

59.1642 

9 

61.1606 

58.3239 

11 

61.8002 

57.6125 

13 

62.4112 

56.9682 

15 

62.9980 

56.3634 

17 

63.5685 

55.7839 

19 

64.1272 

55.2214 

21 

64.6779 

54.6728 

23 

65.2240 

54.1294 

25 

65.7626 

53.5960 

27 

66.2978 

53.0695 

29 

66.8313 

52.5458 

31 

67.3627 

52.0259 

33 

67.8923 

51.5091 

35 

68.4205 

50.9949 

37 

68.9478 

50.4830 

39 

69.4741 

49.9730 

41 

7C.0000 

49.4C4S 

43 

70.5249 

48.9582 

45 

71.0497 

48.4530 

2  _  2(N2  +  N  •*  1 )  (2N  ♦  1) 

^NO  N(N  +  1)  ’  1  ' 

and 

EN/k  =  2.06844N  (N  +  1).  (7) 

The  complete  expression  for  the  absorption  coefficient  is 

a(f,  p,  T)  «  CpT-3  f2  T  An  (8) 

N 

where  d  is  the  atmospheric  pressure  and  T  is  the  absolute  temperature.  For  f  and  Af  in 
gigahertz,  p  in  millibars,  T  in  degrees  Kelvin,  and  a  in  decibels  pir  bolometer,  C  =  2.0058. 
For  a  in  decibels  per  nautical  mile,  C  =  3.7148.  (One  nautical  mile  is  1852  meters 
exactly.) 

Van  Vleck  (9)  estimated  that  Af  was  approximately  0.6  GHz  at  1  atmosphere  pressuri 
and  standard  sea-level  temperature  and  that  it  was  proportional  to  pressure  and  inver-  >ly 
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proportional  to  the  square  root  of  absolute  temperature.  Later  workers  have  suggested, 
based  on  experimental  evidence,  that  the  dependence  on  pressure  and  temperature  is  more 
complicated.  The  Rebel  Mitchell -Carter  (12)  modification  of  the  Mccks-Lilley  (11)  model 
for  this  dependence  is  given  in  terms  of  altitude  h  above  sea  level: 

Af  =  g(h)  ( p/p0 )  (T0/T),  (9) 

where  p0  =  1013.25  millibars  (760  torr),  T0  =  300°K,  and 

g(h)  =  0.640,  0  <  h  <  8  km,  (10) 

g(h)  =  0.640  +  0.04218(h  -  3),  8  <  h  <  25,  (11) 

g(h)  *  1.357,  h  >  25  km.  (12) 

These  are  the  equations  and  constants  that  have  been  coded  in  Fortran  in  Subroutine 

ALPHA,  listed  in  the  last  section  of  this  report,  for  computation  of  ALPH02,  ihe  oxygen 

absorption  coefficient. 


ABSORPTION  COEFFICIENT  FOR  WATER  VAPOR 

The  absorption  due  to  water  vapor  in  the  frequency  region  below  100  GH can  be 
separated  into  two  components,  one  due  to  the  water-vapor  resonance  at  22.2?')  GHz 
and  the  other  due  to  the  residual  effect,  of  the  many  resonance  lines  above  100  GHz.  The 
absorption  coefficient  due  to  the  22.235-GHz  line,  which  will  be  denoted  o22>  ,ias  ■)oen 
computed  from  equations  given  by  Liebe  (14).  The  following  equation  is  a  combination 
of  his  equations  2a,  6,  7,  and  8: 

7  /  2 

*22  =  (~)  (1° 'ogio  e)  S°  PW  (-^)  e2.l 44(1  -3°° /T )  p  dB/km,  (13) 

where  the  symbols  are  defined  as 

c  =  2.998  X  10s  km/sec, 

S°  =  13.92  Hz/torr,* 

pw  —  partial  pressure  of  water  vapor  (torr), 

T  —  temperature  (degrees  Kelvin), 

F  —  line-shape  factor. 

As  was  done  for  oxygen,  F  is  taken  to  he  the  Van-Vleck-Weisskopf  function  (here  includ¬ 
ing  the  factor  f/fr): 


•The  theoretical  value  (14)  is  14.33  Hz/torr,  hut  the  author  has  been  advised  by  Dr.  Liebe  (private 
communication  dated  June  ii,  1972)  that  this  value  involves  approximations  and  that  the  experimentally 
determined  value  13.92  is  probably  more  accurate 
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F  - 


_Af _ 

f)2  +  (Af)2  +  (Tr 


Af 

f)2  +  (Af)2 


whore  fr  =  22.235  GHz. 


(14) 


The  equation  given  by  Liebe  for  Af  is 

Af  =  17.99  X  10“3  (pw  [(300/T)  +  0.20846(pt  -  pw)j  (300/T)"-63].  (15) 

where  again  pw  is  the  water -vapor  partial  pressure  in  torr  and  pt  is  the  total  atmospheric 
pressure  in  torr;  T  is,  as  before,  the  temperature  in  degrees  Kelvin.  The  atmospheric  pressure 
p  in  millibars  is  given  by 


p=  1.33322  pt,  (16) 

and  the  water-vapor  partial  pressure  in  torr  is  given  in  terms  of  the  water-vapor  density  p 
in  grams  per  cubic  meter  (for  atmospheric  concentrations)  by 

pw  -pT/288.75.  (17) 

Combining  all  of  the  numerical  constants,  with  f  expressed  in  gigahertz,  p  in  m;'!ibars, 
end  pw  in  torr,  results  in 

a  =  2.534  X  1CT3  lfpw  (300/T)" ' 2  e2  1J4,;  ‘  nuu'*  F]  dB/l;m  (IS) 

and 

Af  =  17.99  X  lO-3  {pw((300/T>  +  0.20846  (0.75p  -  nw))  (300/T)°  63 }.  (19) 

These  are  the  equations  and  constants  that  are  Fortran  coded  in  Subroutine  ALPHA 
for  the  computation  of  ALPH22,  the  absorption  coefficient  for  the  22.235-GIIz  water- 
vapor  line,  except  that  the  constant  2.534  in  Eq.  (18)  has  been  increased  to  4.693  to 
give  the  absorption  in  decibels  per  nautical  mile.  The  water-vapor  parameter  actually 
employed  in  the  model  atmosphere  is  the  water-vapor  density  p  (grams  per  cubic  meter) 
rather  than  the  partftd  pressure  in  torr;  it  is  converted  to  partial  pressure  by  Rq.  (17). 

The  residual  effect  of  water-vapor  absorption  lines  above  100  (1Hz  has  been  assumed, 
following  Ulaby  and  Straiton  (19)  and  others:* 

arcS  =  CR  (P/P0)  (P'Pot  <T0/T)r'/2  ( f / f 0  ,  (20) 

where  p  is  the  water-vapor  density  (grams  per  cubic  meter),  p  is  the  atmospheric  pressure, 
T  is  the  temperature,  p<),  p0 ,  and  Tq  are  the  zero-altitude  values  of  p,  p.  and  '1',  and  f0 
is  a  reference  value  of  f.  The  constant  Cr  can  he  taken  from  experimental  results  at  any 


•Ulaby  and  Slraiton  aUribun-  Ibis  exponent  for  llic  dependence  oil  temperature  to  Van  V I  e  c  ’< .  although 
it  is  not  given  cxplicity  in  Van  Vleck's  paper  on  water  vapor  absorption  (10).  It  is  probuh'y  deduced 
from  the  footnote  on  page  -120  of  bis  paper  on  oxygen  absorption  (0) 
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suitable  frequency  f0.  Based  on  measurements  reported  by  Straiton  and  Tolbert  (13), 
the  value  Cr  =  0.4  dB/km  at  f0  ~  100  GHz  has  been  adopted  here.  For  p  in  millibars, 
taking  p0  =  1023.23  millibars*,  T0  =  288.1 6°K,  p0  =  7.5  grams  per  cubic  meter  (the 
value  for  the  standard  atmosphere  which  will  be  described  in  a  later  section),  and  fg  = 

100  GHz,  this  formula  becomes 

o rt.s  =  7.347  X  10  3  ppT-5/2f2  dB/km.  (21) 

This  is  the  expression  that  has  been  coded  in  Fortran  for  the  absorption  coefficient  (named 
ALFRES  in  Subroutine  ALPHA)  due  to  the  residual  effect  of  the  water-vapor  lines  above 
100  GHz  except  that  the  constant  7.347  X  10“ 3  has  been  changed  to  1.361  X  10~2  to 
give  the  abjorption  in  decibels  per  nautical  mile. 


DISCUSSION  OF  ABSORPTION-COEFFICIENT 
FORMULAS 

The  total  absorption  coefficient  of  the  troposphere  is  the  sum  of  the  oxygen  and 
water-vapor  coefficients,  and  the  water-vapor  coefficient  is  in  turn  the  sum  of  the  coefficients 
for  the  22.235-GHz  line  and  the  residual  effect  of  lines  above  100  GHz.  Of  these 
components  of  the  total  absorption,  the  one  that  has  the  poorest  theoretical  basis  is  the 
residual  effect  of  the  above-100-GHz  (millimeter-wave  and  infrared)  water-vapor  resonances. 
As  has  been  well  known  for  some  time  ,'10,14,18)  the  application  of  the  theoretical  line 
strengths  and  the  Van-Vleck-Weisskopf  collision-broadening  formula,  Eq.  (14),  give  absorp¬ 
tion  values  that  are  too  low  by  a  factor  of  at  least  5.  Other  line-shape  formulas  have 
neen  proposed,  but  no  fully  satisfactory  theoretical  solution  of  the  problem  has  been 
advanced,  to  the  author’s  knowledge.  Consequently  the  empirical  formula  and  constant 
used  here  are  believed  to  be  as  good  a  solution  as  car.  be  applied  at  present. 

The  use  of  the  Van-Vieck-Weisskopf  line-shape  formula  (20)  for  oxygen  and  for  the 
22.23F  GHz  lines  is  similarly  suspect  for  absorption  prediction  at  frequencies  far  removed 
from  the  resonance  frequencies.  However  this  matter  is  not  so  serious  for  these  lines, 
because  at  frequencies  far  removed  from  the  22.235-GHz  line  its  contribution  is  dwarfed 
by  the  contributions  of  the  oxygen  and  higher-frequency  water-vapor  absorption  lines 
and  hecause  the  oxygen-line  contributions  at  frequencies  above  about  90  Ghz  are  similarly 
dwarfed  by  the  contributions  of  the  millimeter-wave  water-vapor  linos,  except  in  a  very 
narrow  region  about  the  isolated  oxygen  line  at  118.75  GHz.  The  oxygen  absorption  is 
also  small  compared  to  the  water-vapor  absorption  in  the  near  vicinity  of  the  22.235-GHz 
water-vapor  line.  At  frequencies  well  below  the  z2.235-GHz  line  the  nonresonant 
component  of  oxygen  absorption  is  the  dominant  factor.  No  criticism  of  the  shape  factor 
for  this  component  (Eq.  (2))  has  been  found  in  the  literature  by  the  author. 

As  mentioned  in  the  Introduction,  the  results  for  oxygen  and  for  water  vapor  are 
presented  in  this  report  separately  as  well  as  in  combination,  so  that  correction  to  the 
latter  can  be  made  for  any  desired  deviation  from  the  standard  7.5-gTam-per-cubic-meter 
water-vapor  content.  It  was  there  stated  that  a  simple  multiplicative  factor  pip 0  could  be 
applied,  where  p  is  the  desired  water  vapor  density  and  “  7.5  is  the  standard  value. 

As  Eq.  (39)  shows,  this  is  not  quite  true;  there  is  a  small  nonlinear  dependence  of  Af  on 


♦This  value  or  p()  is  the  total  pressure,  consisting  of  the  dry*atmospherc  standard  pressure  of  1013.25 
millibars  and  the  wfter-vapor  partial  pressure  of  9.98  millibars  corresponding  to  p  =  7.5  . 
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p  (through  pw)  in  addition  to  tno  linear  explicit  dependence  of  a  on  p,  Eq.  (18).  A 
similar  slight  nonlinearity  would  appear  in  the  residual  contribution  of  the  higher-frequency 
lines  if  an  exact  expression  were  used,  but  again  the  nonlinearity  is  slight,  and  in  view  of 
the  empirical  nature  of  Eq.  (20)  it  was  not  thought  worthwhile  to  include  it. 


MODEL  OF  THE  ATMOSPHERE 

The  absorption  coefficients,  for  which  formulas  have  been  given  in  the  preceding 
sections,  describe  the  rate  at  which  signals  traversing  an  atmospheric  path  will  be  attenuated, 
as  a  function  of  the  pressure,  temperature,  and  water-vapor  density.  To  find  the  altitude 
dependence  of  the  absorption  coefficient,  it  is  necessary  to  have  a  model  of  the  altitude 
dependence  of  pressure,  temperature,  and  water  vapor  in  the  atmosphere. 


Pressure-and-Temperature  Model 

The  model  adopted  here  for  the  pressure  p  and  temperature  T  of  the  standard 
atmosphere  is  the  U.S.  extension  to  the  ICAO  Standard  Atmosphere  (21);  this  is,  for 
the  altitudes  considered  here  (up  to  100,000  feet  or  30.48  kilomeiers),  the  same  as  the 
ARDC  model  atmosphere. 

The  equations  describing  the  temperature  (degrees  Kelvin)  and  pressure  (millibars) 
of  this  atmosphere  in  terms  of  the  geopotential  altitude  hg  (meters)  are* 

T  =  288.16  -  0.0065  h  "i 

>  hK  <  11000,  (22) 

p  =  1013.25  (T/288.16)0  J 


T  =  216.66 

P  =  (226.32/T)  exp  [— ^(hg  -  11000)]  j 


11000  <  hg  <  25000, 


T  =  216.66  +  0.003  (hg  -25000) 
P  =  24.886  (216.66/T)7 


25000  <  hg  <  47000. 


(23) 

(24) 


where 


a  =  5.2561222, 

=  0.034164794 
7  =  11.388265. 


The  relationship  between  geopotential  altitude  hg  and  geometric  altitude  ha  is 


hg  =  rha  /  (r  +  ha), 


(25) 


Reference  21  ,  pp.  6  and  7,  Eqs.  (1  3)  and  (1  !>). 
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where  r  is  the  radius  of  the  earth,  taken  to  '  356,766  meters.  This  relationship,  which 
is  Eq.  (10)  of  Ref.  21,  is  valid  for  the  alti*  nge  considered  here. 

The  absorption  losses  were  computed  in  terms  of  altitude  in  feet  and  deeihels  per 
nautical  mile,  to  conform  to  the  measurement  units  employed  by  Navy  radars.  The  con¬ 
version  from  feet  to  meters  is  1  foot  =  0.3048  meter  exactly;  the  conversion  from  nautical 
miles  is  1  nautical  mile  =  1852  meters  exactly. 

This  model  of  the  atmosphere  is  incorporated  into  the  absorption-coefficient  computer 
subroutine  ALPHA,  as  a  table  of  values  at  75  specified  altitudes,  in  the  form  of  Fortran 
DATA  statements.  The  pressure  is  given  by  the  dimensioned  identifier  PP  and  the  tempera¬ 
ture  by  TT.  The  75  altitudes  (hf)  range  from  0  to  100,000  feet  in  steps  Ah  that  are 
graduated  according  to  the  following  scheme; 

Ah  =  100  ft,  hr  <  2000  ft. 

Ah  =  1000  ft,  2000  <  h(  <  30000  ft. 

Ah  =  2000  ft,  30000  <  hf  <  70000  ft, 

Ah  =  5000  ft,  70000  <  hf  <  10000  ft. 

This  graduation  provides  more  closely  spaced  points  at  low  altitudes,  where  absorption  is 
greater  and  therefore  where  more  accurate  interpolation  of  the  tabulated  values  is  desirable 
The  interpolation  is  provided  by  the  numerical  integration  subroutines  used  in  computing 
the  cumulative  absorption.  As  will  be  discussed  in  detail  later,  the  integration  is  done  by 
a  modification  of  Simpson’s  rule,  so  that  the  interpolation  is  more  accurate  than  would  be 
obtained  by  using  the  trapezoidal  rule  (linear  interpolation). 


Water- Vapor  Model 

No  accepted  standard  profile  of  water  vapor  has  been  universally  adopted.  The  one 
used  here  is  taken  from  a  report  by  Sissenwine  ot  al.  (22).  Their  values  for  the  water- 
vapor  density  up  to  32  km  (104,987  ft)  are  given  in  Table  2.  As  shown,  the  vapor  density 
is  specified  at  2-km  intervals,  and  the  surface  value  is  5.947  g/m3.  The  values  correspond¬ 
ing  to  the  75  altitudes  specified  in  Subroutine  ALPHA  were  found  by  interpolating.  These 
values  were  then  multiplied  by  the  factor  7.5/5,947  =  1.26114,  to  convert  them  to  the 
surface  water-vapor  value  of  7.5  g/m3  that  seems  to  have  been  adopted  as  a  standard  by 
workers  in  the  field  of  water-vapor  absorption.  The  interpolation  was  done  by  use  of  a 
special  computer  subroutine  by  a  method  equivalent  to  using  a  French  curve  on  a  plot 
of  the  data  points.  This  method  is  described  (in  terms  of  curve  plotting)  in  an  NRL 
Memorandum  Report  (23).  The  resulting  values  of  water-vapor  density  are  incorporated 
in  Subroutine  ALPHA  in  the  form  of  a  DATA  statement  for  the  dimensioned  identifier  RR 

Subroutine  ALPHA  also  has  a  COMMON  statement  containing  a  variable  named 
RHOFAC.  This  quantity  is  set  equal  to  1  in  a  DATA  statement,  but  may  be  set  to  other 
values  by  use  of  the  COMMON  statement  in  the  main  program  or  in  another  subroutine. 

All  the  values  of  water-vapor  density  in  the  model  atmosphere  are  multiplied  by  RHOFAC 
in  Subroutine  ALPHA.  Consequently  if  for  example  calculations  of  the  attenuation  and 
noise  temperature  for  a  surface  water-vapor  density  of  15  g/m3  are  wanted,  the  desired 
result  will  be  obtained  by  setting  RHOFAC  =  2. 
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Table  2 


Midlatitude  Mean  Water-Vapor  Densities* 


Allitude  (km) 

Density  (g/m3) 

0 

5.947  X  10° 

2 

2.946  X  100 

4 

1.074  X  100  ! 

6 

3.779  x  10-1 

8 

1.172  X  10-1 

10 

1.834  X  10-2 

12 

3.708  X  10-3 

14 

8.413  X  10-4  | 

16 

6.138  X  10-4  i 

18 

4.449  X  10-4 

20 

4.449  x  10-4 

22 

5.230  X  10-4 

24 

6.138  X  10  4 

26 

7.191  X  10-4 

28 

5.230  x  10-4 

30 

3.778  X  10”4 

32 

2.710  X  10-4 

’From  Ref.  22,  Table  3. 


This  procedure  may  result  in  vapor  densities  that  exceed  the  saturation  values  for 
the  temperatures  assumed.  However  this  is  probably  not  a  serious  violation  of  realism. 

It  means  that  the  temperature  values  should  also  be  adjusted  when  the  water-vapor  density 
is  increased  significantly.  But  since  the  water-vapor  attenuation  depends  somewhat  weakly 
on  the  temperature,  the  inaccuracy  that  results  from  failure  to  readjust  the  tropospheric 
temperatures  to  values  compatible  with  the  increased  vapor  densities  is  probably  slight. 

i 

It  is  noteworthy  that  p  <h)  does  not  decrease  monotonically  in  this  mode!.  It  does 
decrease  up  to  18  km,  remains  constant  to  20  km,  and  then  actually  increases  from  20  to 
26  km;  thereafter  it  again  decreases.  However,  this  nonmonotonic  behavior  is  probably 
not  too  significant  from  the  absorption  viewpoint,  because  it  occurs  at  a  vapor-density' 
level  that  is  too  small  to  contribute  appreciably  to  the  absorption. 


Absorption  Coefficients  in  the  Model  Atmosphere 

The  variation  of  the  absorption  coefficients  at  zero  altitude  in  this  model  atmosphere 
for  oxygen,  the  22-GHz  water-vapor  line,  the  residual  effect  of  higher-frequency  waier- 
vapor  lines,  total  water  vapor,  and  combined  oxygen  and  water-vapor  absorptions  are 
plotted  for  the  frequency  range  100  MHz  to  100  GHz  in  Figs.  1  through  5.  Figure  6 
illustrates  the  behavior  of  the  absorption  coefficients  with  altitude  at  several  frequencies 
in  the  most-used  part  of  the  radar  frequency  spectrum,  from  100  MHz  to  slightly  above 
30  GHz.  These  curves  have  been  machine  plotted  from  computations  made  with  computer 
subroutine  ALPHA,  which  is  listed  in  the  last  section  of  this  report.  The  plots  have  been 
made  in  terms  of  decibels  per  kilometer  to  facilitate  comparison  with  other  similar  curves, 
although  the  curves  of  cumulative  absorption  loss,  shown  later,  are  plotted  to  nautical-mile 
scales. 


FREQUENCY.  MHZ 


Fig.  2  —  Absorption  coefficient  for  the  22.235-GHz-resonance  component  of  water  vapor 
at  standard  sea-level  pressure  and  temperature 
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FREQUtNCf.  MHZ 

Fig.  5  —  Absorption  coefficient  for  oxygen  and  water  vapor  combined  at  standard 

sea-level  temperature  and  pressure 


AlMToDE.  K1 


(c)  oxygen  plus  water  vapor 


I  —  Variation  of  the  total  absorption  coefficient  in  the  standard 
atmosphere  as  a  function  of  altitude  for  several  frequencies 
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EQUATIONS  FOR  CUMULATIVE  ABSORPTION 

To  find  the  total  attenuation  along  a  ray  j>ath  in  the  atmosphere,  it  is  necessary  to 
integrate  the  absorption  coefficient  along  the  path.  This  path  is  formally  defined  by  the 
following  relation  between  range  R  (as  would  be  measured  by  a  radar)  and  height  h: 


R(l»i ) 


(26) 


To  calculate  cumulative  absorption  along  a  ray  path  from  the  earth’s  surface  to  an 
altitude  h,  the  bending  of  the  ray  path  caused  by  refraction  should  be  taken  into  account. 
This  is  a  function  of  the  ray  elevation  angle  as  well  as  the  refractive  index  profile.  The 
refraction  is  significant  only  at  relatively  low  angles.  In  principle  the  refractive  index 
should  be  calculated  from  the  water-vapor/temperature  profile  of  the  atmosphere.  However, 
since  the  refraction  does  not  have  a  major  effect  on  the  absorption,  it  was  decided  to  take 
it  into  account  simply  by  assuming  a  standard  exponential  model  of  the  refractive  index, 
the  so-called  CRPL  Exponential  Atmosphere  (24)  for  surface  refractivity  Ns  =  313, 
given  by 


n(h)  =  1  +  0.000313  e-K\  (27) 

where  n  is  the  refractive  index  and  k  =  0.00004385  for  h  in  feet. 

Taking  the  refraction  into  account,  from  Snell’s  law  the  derivative  in  Eq.  (26)  is 
given  bv 


dR  _  n(h) _ 

dh  /  r  n0~cos  Op  12’ 

V1  n(l  -  h/r0)_ 

in  which  n(h)  is  given  by  Eq.  (27),  n0  =  0.000313,  0o  is  the  initial  elevation  angle  of  the 
ray  (angle  at  h  =  0  relative  to  the  horizontal),  and  r0  is  the  earth’s  radius.  (The  earth’s 
radius  for  the  CRPL  Exponential  Atmosphere  is  taken  to  be  6370  km,  or  2.0899  X  107  ft.) 

The  cumulative  attenuation  along  this  ray  path  is  given  by 

h, 

A(Rj)  =  2  f  a(h)  Jj^dh,  (29) 

*'o 


op  cos  0  o 
n(l  -  h/r0) 


in  which  Rj  =  R(hj),  o(h)  is  the  absorption  coefficient  at  height  h,  and 

ds  =  1  dR 
dh  n  dh 


(30) 


That  is,  ds/dh  is  given  by  Eq.  (28)  with  1  substituted  for  n(h)  in  the  numerator;  R  is  the 
distance  as  measured  by  a  radar,  along  the  ray  path,  and  s  is  the  corresponding  geometric 
distance.  The  factor  2  in  Eq.  (29)  is  used  to  obtain  the  radar  attenuation,  for  two-way 
traversal  of  the  path. 


*••  ■  •  ■*.*/**** 
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The  cumulative  absorption  is  calculated  in  a  Fortran  subroutine  named  ATLOSS,  which 
is  listed  in  the  last  section  of  the  report.  The  derivatives  ds/dh  and  dR/dh  are  computed  in 
tne  short  subroutine  named  DDH;  the  derivatives  are  therein  named  DSDH3  and  DRDH3. 

The  cumulative  integration  is  performed  by  a  modified  Simpson’s  rule  which  is 
described  in  an  NRL  Memorandum  Report  (25).  However  instead  of  the  actual  subroutine 
described  in  that  report  an  equivalent  procedure  is  programmed  in  Subroutine  ATLOSS. 

The  method  allows  the  cumulative  absorption  to  be  found  at  each  of  the  75  ranges  (heights) 
at  which  a  is  specified,  rather  than  only  at  every  other  one  (half  as  many),  as  would 
result  if  the  usual  form  of  Simpson’s  rule  were  employed. 

This  integration  cannot  be  performed  however  at  h  =  0  when  0O  =  0,  because  then 
the  integrand  becomes  infinite.  A  special  technique  is  used  to  approximate  the  integral 
in  a  small  region  about  h  =  0  where  0o  -  0  (actually,  in  this  case,  from  h  =  0  to  h  =  200 
feet);  this  technique  was  described  in  a  paper  by  the  author  (26). 

Subroutine  ALPHA  is  called  by  Subroutine  ATLOSS  for  a  specified  frequency,  and 
the  resulting  75  absorption  coefficients  a(h)  are  used  to  calculate  the  absorption  (using 
Eq.  (29))  at  75  corresponding  ranges  (found  from  Eq.  (26))  along  the  ray  path  for  the 
specified  elevation  angle.  These  points  allow  a  curve  of  attenuation  in  decibels  to  be 
plotted  as  a  function  of  the  range  in  nautical  miles,  with  frequency  and  elevation  angle 
as  parameters.  Such  curves  are  given  for  38  frequencies  in  the  range  100  MHz  to  100  GHz 
for  several  elevation  angles,  as  Figs.  7  through  98.  Below  1500  MHz  only  the  oxygen 
attenuation  is  plotted  (Figs.  7  through  17),  since  water-vapor  attenuation  is  negligible. 

At  1500  MHz  and  above  plots  are  given  for  oxygen  and  water  vapor  separately,  and  for  the 
total  attenuation.  These  curves  can  be  used  for  estimating  tropospheric  attenuation  for 
earth-based  radars,  for  targets  in  the  troposphere.  They  can  also  be  used  for  estimating 
attenuation  for  poi  ,t-to-point  radio  communication  when  one  terminal  is  earth-bused  and 
the  other  is  airborne  (ground-to-air,  ship-to-air,  or  vice  versa),  by  taking  half  of  the  plotted 
decibel  values. 

The  water-vapor  and  total-attenuation  curves  are  plotted  for  the  standard  water-vapor 
density  (at  zero  altitude)  of  ps  =  7.5  g/m3.  The  total  attenuation  for  any  other  value 
of  ps  can  be  obtained  by  multiplying  the  decibel  value  read  from  the  water-vaoor  curve 
by  ps/7.5  and  then  adding  the  result  to  the  decibel  value  read  from  the  corresponding 
oxygen  curve. 

Curves  of  the  absorption  through  the  entire  troposphere,  for  oxygen  and  water  vapor 
separately  and  for  the  total,  have  also  been  plotted,  in  Figs.  99  through  101. 


TROPOSPHERIC  NOISE  TEMPERATURE 

The  troposphere  radiates  thermal  noise,  and  ihc  tropospheric  noise  temperature  as 
seen  by  an  earth-based  radar  or  radio  antenna  is  given  by  the  following  equation  (27): 


oo 

f  O(R)  Tt(R)  ft-  0.2303  /o  a(r)dr  dR, 


=  0.2303 


(31) 
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Fin.  39  —  Oxygen  component  of  the  radar  absorption  loss  for  ray  paths  in  the 
standard  atmosphere  at  7000  MHz 
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Fig.  58  —  Water-vapor  component  of  the  radar  absorption  loss  for  ray  paths 
in  the  standard  atmosphere  at  30  GHz 
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Fig.  61  —  Water-vapor  component  of  the  radar  absorption  loss  for  ray  paths 
in  the  standard  atmosphere  at  32.5  GHz 
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Fig.  62  —  Total  radar  absorption  loss  for  ray  paths  in  the  standard 
at  nosphere  at  32.5  GHz 
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Fig  87  —  Oxygen  component  of  the  radar  absorption  loss  tor  ray  paths  in 
the  standard  atmosphere  at  75  GHz 
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Fig.  88  —  Water-vapor  component  of  the  radar  absorption  loss  for  ray  paths 
in  the  standard  atmosphere  at  75  GHz 
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Fit;.  99  —  Oxygen  component  of  the  earthbascd-radar  absorption  loss  for  a  target  outside  the  troposphere 
(above  100,000  feet  in  altitude)  for  several  elevation  angles  and  the  standard  atmosphere 
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100  —  Water-vapor  component  of  the  earth-based  radar  absorption  loss  for  a  target  outside  the  troposphere  (above 
100,000  feet  in  rltitude)  for  several  elevation  angles  and  the  standard  atmosphere 
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Fi(>.  101  -  Total  earth -based  ■•-adar  absorption  loss  for  a  target  outside  the  troposphere  (above  100,000  feet  in  altitude) 

for  several  elevation  angles  and  the  standard  atmosphere 
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where  dR  is  measured  along  the  ray  path  (Tn,  in  degrees  Kelvin,  is  a  function  of  the  ray 
elevation  angle),  Tt  is  the  thermal  temperature  of  the  troposphere,  defined  by  the  standard 
atmosphere  of  Subroutine  ALPHA,  and  a  is  the  absorption  coefficient.  (The  factors 
0.2303  are  needed  in  this  equation  because  a  is  here  assumed  to  be  in  decibels  per  unit 
distance.) 

The  previously  described  modification  of  Simpson’s  rule  cannot  be  used  for  this 
integration,  because  it  is  an  integration  with  respect  to  range  (R)  rather  than  height  (h). 

The  height  intervals  are  (within  each  of  the  four  height  regions  defined  in  Subroutine 
ALPHA)  uniformly  spaced,  as  required  by  Simpson’s  rule;  but  the  corresponding  range 
intervals  are  not  uniformly  spaced,  because  of  the  nonlinear  relationship  between  range 
and  height  in  the  space  above  a  spherical  earth  and  in  a  refracting  atmosphere  (for  all 
elevation  angles  except  90  degrees).  Therefore  another  special  modification  of  Simpson’s 
rule  was  devised  to  fit  this  situation;  it  is  Fortran-coded  in  Subroutine  Ik'TGRT. 

When  the  absorption  loss  is  gTeat,  as  it  is  at  the  oxygen  resonance  frequency  for 
example,  most  of  the  tropospheric  noise  seen  by  a  ground-based  antenna  comes  from  the 
lowest  layer  of  the  atmosphere;  that  is,  the  integrand  of  Eq.  (31)  is  large  at  R  =  0  and 
decreases  steeply  beyond  R  =  0,  going  quickly  to  virtually  zero.  Because  of  the  extreme 
departure  of  the  integrand  from  a  parabolic  curve  in  this  circumstance,  the  Simpson's 
rule  integration  is  not  accurate.  An  analytic  approximation  to  the  integral  in  this  region 
is  numerically  more  accurate  than  numerical  integration  of  the  “exact”  expression. 

Therefore  such  an  approximation  is  used  for  the  integration  over  the  first  two  height 
(range)  intervals,  corresponding  to  the  first  200  feet  of  atmospheric  altitude.  The  approxima¬ 
tion  is  derived  by  assuming  that  both  a  and  Tt,  Eq.  (31 ),  are  constants  within  those 
intervals.  If  Rj  and  R2  are  the  ranges  at  n  =  100  and  h  =  200  feet,  then  Eq.  (31)  for 
these  intervals  may  be  written 

Ri  _  ,R2 

ATn=Tl2  l  k,e“klRdR  +  Tt  l  k2e~k2R  dR,  (32) 

2  ->o  2  .'R, 

where  kj  =  0.2303a  j ,  k2  =  0.2303a 2,  Tt|  is  the  average  tropospheric  temperature  in  the 
range  interval  0  to  Rj,  Ttz  is  the  average  temperature  in  the  interval  Rj  to  R2,  and  Oj 
and  a2  are  the  corresponding  average  absorption  coefficients.  Integrating  this  expression 
gives 

ATn  =  Tti  (1  -  e_k  1  R>)  +  Tt  (e~k  2  R>  -  e~k2  R*).  (33) 

This  is  the  approximation  employed  in  Subroutine  AT  LOSS  for  the  first  two  range  intervals. 
(The  approximation  is  used  whether  the  attenuation  is  large  or  not,  because  it  is  equally 
accurate  when  the  attenuation  is  not  large.) 

The  noise  temperature  is  calculated  in  Subroutine  ATLOSS,  since  it  is  efficient  to 
calculate  it  concomitantly  with  the  absorption.  The  result  is  reported  as  the  output 
parameter  ATMP.  A  plot  of  the  resulting  values  for  several  elevation  angles  in  the  frequency 
range  100  MHz  to  100  GH/.  is  given  in  Fig.  102.  This  quantity  is  not  computed  and 
plotted  separately  for  the  oxygen  and  water-vapor  contributions,  because  they  are  not 
linearly  additive  as  are  the  absorption  values. 


IITiTil 


Fig.  102  —  Tropospheric  noise  temperature  computed  for  the  standard  atmosphere  at  several  elevation  angles 
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COMPUTER  ROUTINES 

The  Fortran  subroutines  ATLOSS,  ALPHA,  DDH,  and  INTGRT,  which  have  been 
discussed  in  the  preceding  sections  of  this  report,  are  listed  on  the  following  pages.  Their 
calling  sequences,  pertinent  COMMON  blocks,  and  definitions  of  the  parameters  sue  as 
follows- 

SUBROUTINE  ATLOSS  (FMHZ,  ELEV,  ATMP) 

C0MM0N/RGA/RG(75),  ATTN(3,75) 

FMHZ  —  radio  frequency,  megahertz  (input) 

ELEV  —  initial  ray  elevation  angle,  degrees  (input) 

ATMP  —  tropospheric  noise  temperature,  degrees  Kelvin  (output) 

RG  —  75  monotonically  increasing  values  of  range,  nautical  miles,  along 
the  ray  path  from  h  =  0  to  h  =  100,000  ft  (output) 

ATTN  —  75  corresponding  decibel  radar  attenuation  values  (output) 

for  oxygen  (1),  water  vapor  (2),  and  oxygen  plus  water  vapor  (3) 

SUBROUTINE  ALPHA  (FMHZ) 

COMMON/PTR/PP(75),  TT(75),  RR(75),  ALPH(3,75) 

C0MM0N/H20/RH0FAC 

FMHZ  —  frequency,  megahertz  (input) 

PP,  TT,  RR  —  standard  atmosphere  values,  specified  in  DATA  statements 
(available  as  output  via  the  COMMON  block  if  desired) 

RHOFAC  —  A  numerical  factor  by  which  the  water-vapor-density  values  of 
the  standard  atmosphere  are  multiplied.  If  RHOFAC  is  not  set 
to  some  other  value  by  the  user,  the  value  RHOFAC  =  1,  set  by 
a  DATA  statement,  will  apply. 

ALPH(!,J)  —  the  3-by-75  output  array  of  absorption  coefficients.  The  1  subscripts 
have  the  meanings  1  =  1,  oxygen;  I  =  2,  water  vapor;  I  =  3,  oxygen 
plus  water  vapor.  The  subscripts  J  =  1  through  75  correspond  to 
the  75  altitudes  to  which  the  RG  values  of  Subroutine  ATLOSS  also 
correspond. 

SUBROUTINE  DDH(H) 

COMMON/RRG/REF0,  RAD,  GRAD,  U 
COMMON/DRS/DSDH3,  DRDH3,  AN 


70 


L.  V.  BLAKE 


H  —  altitude,  feet  (input) 

REFO,  RAD,  GRAD,  U  —  refractive-index  profile  parameters  transmitted  from 

AT  LOSS 

DSDH3,  DRDH3  —  ds/dh  and  dR/dh  corresponding  to  H  (output) 

AN  —  refractive  index  n(h)  at  H  (output) 

SUBROUTINE  INTGRT  (HI,  H2,  Yl,  Y2,  Y3,  AREA) 

If  XI,  X2.  and  X3  are  successive  points  on  the  x  axis  of  a  cartesian  coordinate 
system  and  Yl,  Y2  and  Y3  are  the  corresponding  y  values,  then 

HI  =  X2  —  XI  and  H2  =  X3  —  X2.  AREA  (output  parameter)  is  the  value  of 
the  integral' 


-X3 

l  y(x)dx,  (34) 

where  y(x)  is  the  second-degree  polynomial  that  passes  through  the  points  Yl, 

Y2,  and  Y3.  It  is  not  necessary  that  HI  =  H2. 

The  subroutine  listings  follow.  (Slightly  different  but  essentially  the  same  subroutines 
are  used  in  a  computer  program  that  calculates  the  maximum  range  of  a  radar;  tins  program 
is  described  Ref.  1.)  The  lengths  of  these  subroutines  (number  of  memory  locations 
required)  are  as  follows: 

Subroutine  Name  Octal  Length  Decimal  Length 

ATLOSS 
ALPHA 
DDH 
INTGRT 


756 

494 

520 

336 

211 

137 

112 

74 

2021 

1041 

Totals 
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Su'MSuT  I  Nfc  AT1_0SS<F'1HZ.«-LfcV,AT‘ip) 

Cii'iHe^  /PTR/  HP(  75),  17(75),  HR < 75 ) , AtPH<3, 75 ) 

Cl-HhON/RGA/  Ru(  76  )  ,  AT  !  M(  3,  75  ) 

DiMfcMS  t  f,N  166(4).  ntl  M(4  ) 

CoHHflN  /RHj/RtPQ.WAJ.ljRAP.J 

Ci^  1 G  N  /  P  R  5  /  TbHHJ , 6RPH7 , AN 

PArA(^G(l)*0,). (ATT. ^(^,1)*^,), (ATTN (?,1)sO.)«(*7TN(3i1)*6.) 
Data  <RFF£)*,0(j03lJ>,<KAP  =  203RB«5').l3>.(GRAr!s,o0004385) 

ruTA<HAST*0.)»<£LASTsiu0.>i<CONST».?A(52585> 

data  (166*10, 14,  10. 3).  (liPLR*l0O..  1000.  .2000.  ,  50Q0 , ) 

A!Tl(YV)=MC?»(1.25»Vl*2.*V2-.?5*YV) 

AtT2<YY)sfAC2«(-.25»Yl*/.»Y2»lt25*YY) 

Rul(0H)s)ACl«(l.?3*l)RUWl»<;.*,)«plH2'.25*nR) 

RO?(DK):fACl»(-.2>»DB|JM]*2.«0HTH?*l.?b«DP) 

I F  (fMH2  .  to  .  TlAST  .  ANT  ,fLt  V  .fcO.f  LAST  )  riFTU»N 
FlAST  *  tuFV 
1  O  1  S  1 

TnFTAst'_FV/67. 295/796 
S'«:SIM  iHfcTA) 

CbsC06( THfc  T A  ) 

S  5  *  S  N  •  5  N 
RR 1 s 1 , *RF  F  0 

U* ( RPi*SN  )  •  -  Rfc  F  0  *  RFF  0 

IF  < f M H 2 . fcO ,  F  L AS1  1  GO  T ft  55 

CALL  AlPKA(FmhZ) 

Fl  AST*f'MM2 

6  M:0, 

R.nGsO  , 

A  !  T6M*ATTEN?  =  0  , 

K*-l 

MM  I  ^  *  0  , 

I  F  (  6lE  V  .  c2, 0  .  )  M«I  N«1  .6 -9 
CALL  DDH(HM]N> 

OKPHlxPROHj 
05OHl*nS0M.3 
A  N  1  s  A  N 

TKl.ALPM(3,l)»tT(l> 

TtMP  «  0 , 

Yls ALPH ( 1 , 1 ) • JSD>U 
Y 1  1  *  Y  1 

Y1?*ALPH(2 , 1 >»DSDhl 
Po  60  J  =  1  ,  4 

FaC1ju6lH(J)/(3. *6076. 1155) 

F AC2*2 , •' ACl 
I  MA  <:  I  66 ( J ) 

Dli  61  III,  IMAX 
K  =  *«? 

HS  H  ♦Del.  (  w  > 

HI  *H 

C  A  i  L  UOH(H) 

n-inH2*nRDh3 

D3PH2aDsnH3 

At?sAN 

HsO tOtl H [ j ) 

Call  I)CH(h) 

Yi s  A  l  PH ( 1 , K ♦ 1  ) •DSDH2 
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<21  =  Y2 

Y22=AlPH(2,K*1 
y3:Alph< 1,K*2)»DSCHJ 
Y-IPsALPh^.K^^DSOHS 
A  N  5  *  A  i* 

!►  (ElEv  ,LT.  1.  ,ANO.  H  ,  i_T  ,  201. >  5,  6 
0  Cl»CS»CS*(1./RaD-HEF(i»RRAU/RP1) 

CU1  *  1  .  /  (  CC»60  >6 , 11*5&  ) 

Pk90s2 . .CC*Hj 

C  FuLLG* 1 NG  IS  APPLICATION  WEQ'lIPED  NEAR  TWETA*0  *ND  HaO  F®ft  RANGE 
C  CALCULATION,  range  lb  CALCULATED  thus  FOR  h  s  100  AND  Ht  200  whEn 
C  cuPvAlIfiN  ANC.lE  is  lEsR  than  1  DEGREF. 

Ri\r,sCCl*PR0O/(SQRTE  (PkOd.5S)*SN) 

OS1 *  RNG 

r  approximate  attenuation  is  range  (Two-way)  times  average  value  e® 

C  ALPHA  In  THE  range  i  n  I  e  h  V  *L  , 

aITEN1*RNG»(ALPH(1.1)  ♦  ALPHtl.2)) 

AT  TgNisRNGM  ALPH(2. 1  )  ♦  ALPH  (  2  «  2  )  ) 

C  radar  Range  IS  GEgmET h I C  RANGE  TIMES  AVERAGE  REFRACTIVE  INDEX, 

RnO  =  R,*G »  (  HPl  .  1.  *  PeFG»fcXP< -GRaD»H1  )  )  •  ,5 
I  S  I  G  s  2 
GO  TO  7 

6  DsrRGl(DRCH3) 

RnGsRNG»DS 

Db 1 =  US/ ( ( ANl »AN?) • ,5) 

ATTgNl*ATTfNl*ATTKY3) 

Y  1  -  Y  1  2 
Y2l Y?2 

A ! TEN2*ATTEN?  ♦  A|Ti<r32) 

7  RG(2*1):RNG 
a:tn<i,k»i)*  ATTENl 
Al TN(2,K*l)sATTEN2 

A  1  T,N(  J  .  K*1  1  *AT  TENl  «  A  T  1  E  N? 

GO  TO  (10,11)  IS1G 
U  PkOD=2,«CC«H 

RNr,sci;i»PR0O/(SORrf  (proi.*ss)«sn) 

DbAsRNG- DS1 

ATT6NX*RNG»(AlPH(1,1)  ♦  ALPH(l.J)) 

AT T6N2rRNG« ( ALPH(2 . 1 )  *  ALPH(2,3)) 

RNGiRNG* ( RP1  ♦  1.  *  PtPo»fcXP(*CRAD«H  ))  *  ,5 
1  S!S  =  1 
Gu  TO  12 

10  DssRG2(0R0H3) 

RN  ",  s  R  N  G  ♦  0  S 

Db5srs/(  ( AN2«AN3  )• .> ) 

YlrYl 1 

Y  2  s  Y  2 1 

A  1  T  E  N 1  *  ATTENl  ♦  ATT2(Y3) 

YlrY12 

Y  2  s  Y  22 

AITEN2=  ATTEN2  •  AT  T? ( YO?  ) 

1?  Rj(K*2)-RNG 

A ' Tn ( 1 ,*  *2  )  s A  I  TEN1 
a;  Ts(2.K42):AITEN2 
ATTN(J.k*2)5AITEN1  ♦  AT  TEN? 

Ai_RSSsii-.»»(-ATTN(3,K»?)»,n>) 


non 
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TP3aAlPh(,J.>W>»TT<iv*2>»AL':'SS 

Cil  *  <RG(K*l)-sri(K))/((ANl*AN7).,rj) 

Db?  S  (  HG(K»2  ) -RG  (  K*1  )  >/(  (  AN2*AN3  )  •  .  *>  ) 

1 ►  (K.EO.l)  70.71 

APPROX  1  MA  T  ION  EHPtl'YtLl  IN  PLACE  Of  FIHST  l  N  T  u  «  A  7  1  *  N  SY£D, 

IMcD  Tg  GIVE  V  A  l  l  n  Re«uiTS  IN  Hi  GH-A  I  TEND*  T  ldN  r  a  R  t  ?  •  analytic 
A  P  P  ROX I MAT I0N  STARTS  AT  STATEMENT  70. 

/0  CtX*0,5»CQNST««ALPH(3.1).AuPH{T,2)) 
CkVtOt9*CeNS>'*(ALPH(S»?)»ALPH(^.3n 
AuOSl «fc  XPF  (*ChX*0Sl  ) 

AtPS2  =  ExPF  < -CtY.nsi  ) 

AL0S3*ExPf  (•CfcT«(l,'Si*i)S2)  ) 

D!EmP*(o.5/C0NST>»((TI  (1)*TT<2))*»1 , • A  L  0  S 1  1  ♦  (lT(P)*TT(3))*fALi)S2 
1  -AUPS3)) 

GO  T  0  72 

/I  Al0$S«1 0 . • • ( - aT TN < 3  ,  y ♦ 1  )  • . P5  ) 

Ti-2s4LPH<3,K*l>»TT<X*l5«ALPSS 

Call  lNTGRT(Dsi,0S?,TPl.TP2,TP3,nTEMP> 

/ 2  T  t  HP  s  TFMP»DIEHP 
0RPhi*0wPh3 
T  P  1  *  T  p  3 

YlsYU«Y3 
Y12SY32 
AM  a  AN3 
01  CONTINUE 
6  P  CONTINUE 

Al mPiTEmp»CGnsT 
EnH 
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C 


SlRHGcTlNfe  AL?6A<8  K -< 2  ) 

C  E  M  M  f  N  /0T9/  f •P(7a),  TMJili  hp(7i  )  ,  AL^Hf  J>  7*> ) 
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